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The renormalization group flow of an integrable two di- 
mensional quantum field theory which contains unstable par- 
ticles is investigated. The analysis is carried out for the Vira- 
soro central charge and the conformal dimensions as a func- 
tion of the renormalization group flow parameter. This allows 
to identify the corresponding conformal field theories together 
with their operator content when the unstable particles van- 
ish from the particle spectrum. The specific model considered 
is the S'[/(3)2-homogeneous Sine-Gordon model. 
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The study of two-dimensional quantum field theories 
(2D-QFT) has turned out to be a fruitful venture since 
almost three decades. In particular when exploiting inte- 
grability many non-perturbative methods have been de- 
veloped over the years. Besides the challenge to under- 
stand the underlying mathematical structures and the 
intriguing physical applications in two dimensions itself, 
e.g. to describe measurable quantities of carbon nan- 
otubes the ultimate goal is to extrapolate ones find- 
ings to higher dimensions. In particular for the cele- 
brated c-theorem of Zamolodchikov , which originally 
describes the renormalization group trajectory of a func- 
tion which at the renormalization group fixed point cor- 
responds to the Virasoro central charge, various coun- 
terparts have been developed in higher dimensions, e.g. 

Fairly recently a class of massive integrable quan- 
tum field theories, the homogeneous Sine-Gordon models 
(HSG) [Q, has been proposed introducing the feature of 
possessing unstable particles inside its particle spectrum. 
Despite the fact that theories containing resonances have 
been treated before in the context of two-dimensional 
massive quantum field theories, e.g. the HSG-models 
are somewhat special since they constitute the first ex- 
amples of theories which admit a well-defined Lagrangian 
description. In general the HSG models are associated to 
integrable perturbations of G-parafermions of level k [|j , 
i.e. WZNW-coset theories of the form Gk/U{lY with £ 
being the rank of a compact Lie group G. As free param- 
eters the model contains £ different mass scales and i—1 
different scales for the resonance parameter a which en- 
ter the Breit-Wigner formula [Q. In general an unstable 
particle of type c is described by complexifying the phys- 
ical mass of a stable particle by adding a decay width Fg, 
such that it corresponds to a pole in the S-matrix as a 
function Mandelstam s at s = = (Mg — iTc/2)'^ (for 



a more detailed discussion see e.g. [||). As mentioned in 
whenever Mg ^ Fg, the quantity A/g admits a clear 
cut interpretation as the physical mass. However, since 
this assumption is only required for interpretational rea- 
sons we will not rely on it. Transforming as usual in 
this context from s to the rapidity plane and describ- 
ing the scattering of two stable particles of type a and b 
with masses rua and rrib by an S-matrix SatiO) as func- 
tion of the rapidity 6, the resonance pole is situated at 
Oft = a — ia. Identifying the real and imaginary parts of 
the pole then yields 
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m'^^ + ml + 2m am h cosh a COS a (1) 
2m,amb sinh \a\ sin a . (2) 



Eliminating the decay width from (|l|) and (g), we can ex- 
press the mass of the unstable particles Mg in the model 
as a function of the masses of the stable particles TOq, mh 
and the resonance parameter a. Assuming a to be large 
this gives 
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One recognizes the occurrence of the variable me'"'!/^, 
which was introduced originally in ^ in order to describe 
massless particles, i.e. one may perform safely the limit 
m — > 0, cr — > oo and one rnight therefore be tempted to 
describe flows related to (^ as massless flows. In jl^] 
the relative mass scales between the unstable and stable 
particles and the stable particles themselves were investi- 
gated by computing the finite size scaling function from 
the thermodynamic Bcthe ansatz (TBA). A consistent 
physical picture was obtained for the overall identifica- 
tion of the flow between different coset models. It re- 
mained, however, an open question how to identify the 
operator content. In general this question is left unan- 
swered in the context of the TBA. For theories with cer- 
tain properties, it is sometimes possible to determine at 
least the dimension of the perturbing operators by in- 
vestigating periodicities in the so-called Y-systems [ pT] . 
Resorting to a different method, namely by appealing to 
sum rules which are expressible in terms of correlation 
functions, the major part of the operator content was 
successfully identifled for some of the HSG models . 
The purpose of this manuscript is on one hand to confirm 
and refine the TBA results by the latter method, i.e. by 
investigating the renormalization group flow described by 
the Zamolodchikov c- function . We will precisely study 
the onset of the mass scale of the unstable particles and 



1 



investigate how a particular coset flows to another one. 
On the other hand, we study in addition the flow of the 
operator content of one conformal field theory to another 
one by exploiting the flow provided by the A-sum rule of 
Delfino, Simonctti and Cardy [ p^ . 

Denoting by r the radial distance and by i = In the 
renormalization group parameter, the functions c{t) and 
A(t) were defined in [g| and ||l3|, respectively, obeying 
the differential equations 



dt 
dA{t) 

dt 
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The r.h.s. of these equations involve the two-point cor- 
relation functions of the trace of the energy-momentum 
tensor and an operator O, which is a primary field 
in the sense of ||l^. In general these equations are in- 
tegrated from t = — ootoi = oo and one consequently 
compares the difference between the ultraviolet and the 
infrared fixed points. In order to exhibit the quantitative 
onset of the mass scale of the unstable particles we inte- 
grate these equations instead from some finite value to to 
infinity. Restricting our attention to purely massive the- 
ories we use the fact that for those theories the infrared 
central charges are zero, such that 



drr^ (e(r)e(O)) 



(6) 



Instead of the integral representation (||), the c- function 
is equivalently expressible in terms of a sum of correla- 
tors involving also other components of the energy mo- 
mentum tensor [|j . In deriving (^ these terms have been 
eliminate by means of the conservation law of the energy 
momentum tensor. We find (^) most convenient. The 
flow of c(ro) will surpass various steps: Starting with 
tq = the theory will leave its ultraviolet fixed point and 
at a certain definite value, say tq = r^, the unstable par- 
ticle will become massive such that c(ro > r^) can be as- 
sociated to a different conformal field theory. It appears 
natural to identify the mass Mc as the point at which 
c(ro) is half the difference between the two coset values 
of c. As a consequence of @j we may relate the masses of 
the unstable particles at different values of the resonance 
parameter a, a' and expect Mc{ru,a) = Mc{r'^,(j'). We 
will employ the latter equality evaluated in the form (|^) 
not only as a consistency requirement, but also as a con- 
firmation of the fact that the renormalization group flow 
is indeed achieved by to — > tq to. Increasing vq further, 
the energy scale of the stable particles will eventually be 
reached at, say at rg = r^, rf,, . . . , r„. Depending on the 
relative mass scales between the stable particles these 
points may coincide. Finally the flow will reach its in- 
frared flxed point c(ro = Vir) = 0. 

Likewise we can integrate equation (|^) 



A(ro) 
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which allows to keep track of the manner the opera- 
tor contents of the various conformal fleld theories are 
mapped into each other. We used that all conformal 
dimensions vanish in the infrared limit. Fortunately, we 
have (0(r)O(O)) ^ (0(0)) in many applications such that 
the vacuum expectation value (0(0)) cancels often. One 
should note, however, that (j^) is only applicable to those 
operators for which its two-point correlator with the trace 
of the energy momentum tensor is non-vanishing, such 
that one may not be in the position to investigate the 
flow of the entire operator content by means of (|^). 

In order to evaluate (||) and (|^) we have to compute 
the two-point correlation functions in some way. In 2D- 
QFT this is probably most efficiently achieved, by ex- 
panding them in terms of n-particle form factors, i.e. the 
matrix elements of some local operator 0{x) located at 
the origin between a multiparticle in-state and the vac- 
uum denoted by (0|O(0)|F^, (^i)F^, ((?2) ■ ■ • V^M),^ =■ 

F,?'^'-''"(6'i,...,6'„). Here the V^{e) are some ver- 
tex operators representing a particle of species /i. Ab- 
breviating the sum of the on-shell energies as £' = 
S"=i '^Mi cosh6'i, one may write 



oc 

oo „ 

(O(r)O'(0)) = ^ ^ j 



n=l fii...fj,„_ 



ddi . . . dOn _r E 

n!(27r)" ^ 



(8) 



Using this expansion we replace the correlation functions 
in the expression of the c- function c(ro) and the scaled 
conformal dimension A(ro) and perform the r integra- 
tions thereafter. Thus we obtain 
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and 
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We wfll now analyze (|), (0) and (|) for the SU{3)2- 
HSG model. This model contains only two self-conjugate 
solitons which we denote by , "— " and one unsta- 
ble particle, which call u. The corresponding scatter- 
ing matrix was found to be S±± = — 1, S±^{9) = 
± tanh (6* ± cr — i7r/2) /2, which means the resonance 
pole is situated at Or = — in/2. Stable bound states 
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may not be formed. Note that for the corresponding 
value of (7 = 7r/2 and arbitrary a the condition ^ 
is not fulfilled. However, as indicated above this condi- 
tion only helps for a clearer identification of the mass 
parameter. For the HSG-models this condition starts 
to hold when the level is large, which indicates that in 
these type of models this interpretation is in fact a semi- 
classical one. 

A huge class of form factors corresponding to vari- 
ous operators related to this model were constructed in 
p^ , p^ . Labelling an operator by four quantum numbers 
/i, T, t' the general n-particlc solution reads 



0''-",\M+M- 



2S+T,2t + T' 



(01, 



C-' , M+M- 
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2s+T,2t+T' 



2s+t) 



(<^2t+r')" 



"'YlF^'^^ie,,). (11) 



We used here a particular ordering by starting with 2s+t 
particles of the type n = + followed by 2s + t' particles of 
the type ^ = —, collected in the sets M"^ = {±, . . . , ±}. 
Once these expressions are known, all other form factors 
related to it by permutations of the particles may be con- 
structed trivially by exploiting Watson's equations 0, 
see 1 15 1^ for details concerning the HSG-models. The 
functions F^'^^ for all combinations of the fi's are 



F^^{e) = -V2tanh^exp(=F6'/2) 



(12) 



i7r(lTl)±9 G r= dt sin^(('^-^Tg) aV) 
F^^{9)— 2^6 ^ TT ~Jo t sinhtcosht/2 (^l^'j 



with G = 0.91597 . . . being the Catalan constant. The 
{t + s)x{t + s)-matrix 



y^2s+T,2t+T' J ^ 



2(j-i)+^' — — 
^2U-i)+2t+u ' t<l<S + t 



(14) 



has as its entries elementary symmetric polynomials (see 
e.g. ||l^ for properties) depending on different sets of vari- 
ables. We use the notation cr* when they depend on the 
variable x = exp6' associated to the sets and a to 
indicate that all variables are multiplied by a factor ie^"^ . 
The overall constant was computed to 



2S+T,2t + T' 



js(2r+T'+iy+2)2s(2s-2f-T'-l+2r) 



,gS<T(2f+r')/2^ 



r,2t+T' ' 
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where the value of 2t+T' fixed by the lowest non- 
vanishing form factor. In particular we need 



^2^,2t = 0-1 (a^i,. 
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Having assembled all the ingredients we can evaluate the 
expressions (^) and ([l^). We carry out the integrals by 
means of a Monte Carlo computation. For c(ro) we take 
contributions up to the 4-particle form factor into ac- 
count and display our results in figure 1. 
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Figure 1: Renormalization group flow for the Virasoro central 
charge c(ro) for various values of the resonance parameter a. 

Following the renormalization group flow from the ul- 
traviolet to the infrared, figure 1 illustrates the flow from 
the SU{Z)2/U{\f- to the SU{2)2/U{1) ® SU{2)2/U{1)- 
coset when the unstable particle becomes massive. This 
confirms qualitatively the previous observation of the 
TEA analysis Here we also want to compare the 

value of the mass of the unstable particle at different 
points of the resonance parameter a and t^. Taking 
now the mass scales of the stable particles to be the 
same, i.e. ~ m_ — m, we compute the mass of 

the unstable particle according to (^), i.e. Mu{tu,cT) ^ 
TO/\/2exp((|cr| +tu)/2). This means for different values 
of the resonance parameter we may still have the same 
value for the mass of the unstable particle when changing 
tu, indeed we find 

Ms(-30.8, 30) = A/{,(-20.8, 20) = M,i(-10.8, 10). (17) 

Since the flow between the two cosets is smooth and takes 
place over some range of to, we had to select one par- 
ticular point tu- As already indicated in general, it is 
convenient to identify as the point at which c(io) 
is half the difference between the two coset values of c. 
It is clear from figure 1, that since the overall shape of 
the curves between two values of c is identical for dif- 
ferent values of a, any other value in the interval would 
lead to the same results in comparative considerations. 
This also means that when evaluating ( |l7| ) the resulting 
value 0.47m, which apparently violates the energetically 
necessary condition > ma +mi,, should not be taken 
too literally since the point t„ is only chosen because it is 
easy to fix. Equations ( [T7| ) confirm our general assertions 
outlined above. 

For the evaluation of the scaled conformal dimension 
(|lO| ) we proceed similarly. For the solutions correspond- 
ing to the operators Oq q, Oq2 and ©2 0' whose confor- 
mal dimension in the ultraviolet limit was identified (l|] 
to be 1/10, we take up to the 6-particle form factors into 
account. For the former two operators our results are 
presented in figures 2 and 3. 
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Figure 2: Renormalization group flow for the conformal di- 
mension A(ro) of the operator Oq'q for various values of the 
resonance parameter a. 

We observe that the conformal dimension of the op- 
erator Oq'q flows to the value 1/8, which is twice the 
conformal dimension of the disorder operator fi in the 
Ising model. The factor 2 is expected from the mentioned 
coset structure, i.e. we find two copies of SU{2)2/U{1). 
The nature of the operator is also anticipated, since by 

construction F„ of the S'(7(3)2-HSG model co- 

incides precisely with F^^ of the thermally perturbed Ising 
model when one of the sets is empty. It is also clear 
that we could alternatively obtain ( [l7| ) from the analysis 
of A(ro). 



explains the fact that the starting point at the ultraviolet 
fixed point is not quite 0.1. The operators also flow to the 
value 1/8, such that the degeneracy of the 5'{/(3)2-IISG 
model disappears surjectively when the unstable particles 
become massive. 

In comparison with other methods it would be ex- 
tremely desirable to elaborate on the precise relation- 
ship between c(ro) and the finite size scaling function 
of the thermodynamic Bethe ansatz. Also the relation 
to the intriguing proposal in [p^ of a renormalization 
group flow between Virasoro characters remains unclar- 
ified. The analogue of A(ro) still needs to be identified 
in the TBA as well as in the context of . In addition 
one may pose the question whether there exist higher di- 
mensional counterparts of the function A(ro) in analogy 
to the results obtained in Q for c(ro). Concerning the 
specific status of the HSG-models it remains a challenge 
to extend the results to other Lie groups pO| ]. 
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Figure 3: Renormalization group flow for the conformal di- 



mension A(ro) of the operator Og'j 
resonance parameter a. 



for various values of the 



Despite the fact that the explicit expressions for the 



form factors of Oq'I and O2 o differ the values of A(ro) 
are hardly distinguishable and we therefore omit the plots 
for the latter case. We also note the previously observed 
fact [|2j, that the higher particle contributions for the 
latter operators are more important than for Oq'q, which 
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